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ABSTRACT 

/ Cr - ' / ' 

The pressure fluctuations caused by a turbulent 
boundary layer on a solid surface have been obtained 
over the Mach number range 1.33 to 5.00. The root-mean- 
square values of the pressures are found to be proportional 
to the local skin friction. The proportionality constant is 
not a strong function of Mach number in the range investi- 
gated but is larger than the previously obtained subsonic 
values by about a factor of 2. The space — time correla- 
tions with the space separation in the direction of the 
mean flow are characterized by a convection speed U c , 
and this speed falls from 0.8 of the free-steam velocity 
at M = 1.33 to 0.6 at M = 5. The peak value of the correla- - 
tion coefficient falls to one-half for a spacial separation 
of the measuring points of about two-tenths of the boundary 
layer thickness. 


I. INTRODUCTION 

Associated with the irregular motions of a turbulent fluid is a fluctuating pressure field. This 
pressure field is apparent to the unaided observer through both the sound field associated with the turbulence 
and the fluctuating force on a solid surface in contact with the turbulence. Recently there has been an 
increase of interest in the pressures on material surfaces, since when flight vehicles are operated in regimes 
of large dynamic pressures, these pressures can produce significant effects. The random fore es can cause 
fatigue failure in a structure as well as undesirable levels of structural vibration. In addition, these forces 
can produce sound within a structure through the intermediate step of forcing the solid surface into motion. 
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There is, consequently, a definite need to obtain information about the turbulent pressure field, particularly 
for Mach numbers near and above one. With this motivation, a program to obtain measurements of the pressure 
field on a solid surface adjacent to a supersonic boundary layer was undertaken at the Jet Propulsion 
Laboratory, where measurements could be obtained over a large Mach number and Reynolds number range 
in the supersonic wind tunnel. 

Most of the work done previously on this problem, both theoretical and experimental, has been 
limited to the study of the incompressible and the low Mach number cases. For incompressible flows, the 
Navier -Stokes equations give a relation between the pressure gradient and certain combinations of velocity 
derivatives and velocities. The pressure at a point can be computed, therefore, as a line integral of some 
function of the velocity field. A more tractable approach is to take the divergence of the momentum equations 
and obtain a scalar expression for the Laplacian of the pressure. This equation is a Poisson equation, where 
the source terms involve a knowledge of the velocity field; it can be integrated for a particular set of boundary 
conditions by a Green function technique, and the solution is at hand if we can make the correct assumption 
concerning the source terms. Unfortunately, the available experimental results for the turbulent velocity field 
are not sufficient to support a statement about the source terms free from contentious points; hence the 
theoretical results do not contain much more information than could be obtained by a simple dimensional 
argument. Kraichnan (1) and Lilley (2) have studied the theoretical problem. 

Measurements have been made of the wall pressure fluctuations for subsonic turbulent boundary 
layers by Willmarth (3) and Harrison (4) . Their results were used as a guide in setting up these experiments. 

There is a definite advantage in going to compressible flows to get information with which to 
construct a theory, since the ratios of the various derived scales of a shear flow, e.g., the momentum or 
displacement thickness, can be varied independently by varying the Mach number and Reynolds number. Also, 
as the Mach number is increased, some phenomena might be accentuated which could help clarify the physical 
process. 

It would be desirable to obtain measurements of the pressure field throughout the turbulence in order 
to construct a rational theory. Unfortunately, measurements anywhere except at the wall represent an almost 
impossible problem since any pressure-sensitive device in the flow responds both to static pressure and 
local velocity fluctuations. 
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II. AERODYNAMIC EQUIPMENT 


The measured flow fields were produced in the JPL 20-in. supersonic wind tunnel. This continuous 
flow tunnel has a flexible plate nozzle which permits Mach numbers up to 5.0 to be attained. The top 
and bottom surfaces of the 18 x 20-in. cross-section working area are made up of flexible plates, whereas 
the side walls are plane upstream of the throat through the working section. The nozzle produces a uniform 
flow in the working section to within ± 1% of M over the whole Mach number range. The settling chamber 
contains an array of screens to act as turbulence reducers, and measurements have shown that the turbulence 
in the main stream of the working section is almost wholly caused by sound radiated from the turbulent 
boundary layers on the tunnel wallsf^ which is the irreducible minimum unless these boundary layers can be 
removed. 

The specific flow fields investigated were the boundary layers on a flat plate mounted in the tunnel 
and the boundary layer on the flat side wall of the working section. The mean properties of these boundary 
layers closely approximated the properties of equilibrium, adiabatic, flat-plate boundary layers; therefore, 
the boundary layers are assumed to be uniquely characterized by the free-stream Mach number and by a 
Reynolds number based on some thickness. 

The air supply of the tunnel was dried sufficiently to prevent water condensation. The tunnel 
stagnation temperature was approximately 100°F; stagnation pressure ranged from about 20 to 300 cm Hg. 

The flat plate used for these tests was 1 in. thick, 18 in. wide, and 33 in. long. It spanned the 
tunnel along the horizontal midplane of the working section. The plate was flat on the surface upon which 
the measurements were taken, the opposite surface being beveled at 24 deg at the upstream end to form a 
sharp leading edge. The leading edge was perpendicular to the stream direction with a radius of curvature 
of a few thousandths of an inch. A row of small holes through which air could be blown to effect the transi- 
tion of the boundary layer was located near the leading edge of the test surface. The plate installed in the 
tunnel is shown in Fig. 1. 

The major problem in measuring the pressure fluctuations is the isolation of the pressure pickup 
from the accelerations associated with the vibration of the tunnel. Because of variations in tunnels, flow 
fields, and transducers, there is no general solution for this problem. In these tests it was found sufficient 
to support the flat plate in the tunnel with four bolts extending through the side walls in rubber bushings. 
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This mounting scheme together with a suitable transducer, described in Section III, resulted in an output 
due to pressure significantly larger than that from acceleration for a substantial range of flow conditions. 

In order to obtain measurements of the pressure fluctuations on the side wall, it was found necessary 
to mount the transducer on a disk that could be inserted into a hole in the side wall with its surface flush 
with the surface of the side wall (Fig. 2). The disk had to be supported at another point in the tunnel 
structure where vibrations were at a minimum. The area around the disk was sealed with rubber. 

The flat plate configuration was used in the early stages of these experiments and had the advantage 
that both a laminar and a turbulent boundary layer could be produced over the transducer for the same tunnel 
conditions. Also, a large variation in the boundary layer Reynolds number at the transducer site could be 
obtained by varying both the tunnel pressure and the transducer location. The size and frequency requirements 
for quantitative measurements on the plate were so extreme, however, that as soon as the side wall mounting 
problem was solved, most of the data were obtained there. The size and frequency requirements were eased 
by a factor of about five compared to the plate. 
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III THE PRESSURE TRANSDUCER 

The three major requirements for a pressure transducer are that it have a sufficiently large frequency 
range, that it be sensitive enough so that the output is larger than the electrical noise of the associated 
equipment, and that it be small enough to approximate a point compared to the relevant scales of the 
pressure field. 

It is known that the spectrum of the velocity fluctuations in a supersonic turbulent boundary layer 
shows significant energy out to a frequency range of about 5 U^/8, where U ^ is the free-stream velocity 
and 8 is the geometrical boundary-layer thickness — the distance from the wall at which the mean velocity 
in the boundary layer becomes equal to the free-stream velocity. With the stagnation temperature that was 
used, the free-stream velocity was of the order of 2000 ft/sec. On the plate, 8 ranged from 0.1 to 0.4 in.; on 
the side wall, it varied from 1 to 2 in. Therefore, assuming the pressure scales were comparable to the 
velocity scales, the response of a transducer has to extend to 600 kc on the plate and to 150 kc on the 
sidewall. The size of the transducer should be smaller than 0.2 8 or ideally smaller than 0.02 in. on the 
plate. The transducer (Fig. 3, 4) developed for use within these limits used a small lead zirconate element^ 
as the electromechanical transducer. This particular configuration is the result of long development and 
yielded meaningful pressure measurements for the particular set of conditions encountered. 

The design features of the particular mounting are as follows: 

1. The surface of the pressure-sensitive element is flush with the wall; there are no open 
spaces on the surface. It was found that leaving an air gap (of approximately 0.002 in.) 
around the element in order to isolate it from vibration resulted in a strong extraneous 
output whose frequency corresponded to that of a Helmholtz resonator of the volume and 
area present. This frequency fell within the range of interest, i.e., less than 600 kc, and 
the air gap was eliminated. 

2. Only a single piece of PZT-5 was used. All attempts at making a layered structure, 
either to increase the sensitivity or to cancel some accelerations, resulted in extraneous 
signals in the frequency range of interest (1 to 600 kc). The material used was either a 
cut-down piece of PZT-5 disk, or, for thin pieces (0.020-in. thick), a piece of the PZT-5 
element of a Clevite TF-01A electromechanical transformer. PZT-5 has one of the 

1 PZT-5, manufactured by the Clevite Corp., Cleveland, Ohio. 
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highest Curie points of any commercially available piezoelectric ceramics, so that the 
leads could be soldered directly to the silver layer on each face. 

3. The major resonant frequency of the transducer is determined by the thickness of the 
piezoelectric material, and can be increased to well above 600 kc by using material of 
sufficiently small thickness (0.02 in.). 

4. The physical area in contact with the flow can be reduced to a circle 0.015 in. in diameter. 
Transducers with sensitive elements of various sizes were produced, with thicknesses 
ranging from 0.02 to 0.1 in. and diameters ranging from 0.015 to 0.1 in. 

5. Any motion of the leads from the piezoelectric element can produce a spurious output 
signal if the leads are near a conducting surface (capacity change). These leads were 
imbedded in epoxy resin and held rigidly until they could be brought away from the 
structural surface. 

6. In order to get sufficient isolation of the transducer from the vibration of the plate or wall 
disk, it was mounted in rubber. The vibrations of the transducer on this support gave the 
low-frequency limit to the measurements, and in order to make this frequency as low as 
possible, a lead cylinder was used to back up the ceramic element. The size of this 
cylinder was controlled by the available clearance space. 

When the transducer is considered as a unit, the following remarks can be made. The stresses in the 
sensitive element that are interpreted as a pressure signal on the exposed surface can actually come also 
from pressures exerted on the epoxy surface and transmitted through the lead or rubber. That this energy was 
small compared to the direct action of the pressure on the ceramic element was shown by transducers in which 
the epoxy resin had been hollowed out. Only the epoxy surface and edge were left; consequently, the lead 
backing was somewhat decoupled from the surface. In test situations, the outputs of these transducers were 
no different from outputs of transducers with solid epoxy buildups. The possibility that an intense standing 
wave pattern existed in the backing block was also checked by using different backing materials (glass and 
steel) and by modifying the shape of the backing block (grooving). No effect was observed with any of these 
modifications. For some reason no measurement showed a frequency that could be attributed to acoustic waves 
in the backing block. The strongest vibration associated with the transducer itself was the low-frequency 
resonance of the lead backing on its rubber support, which occurred at around 500 cps. Of course, a particular 
transducer might show unexplained signals at any frequency, but these could be attributed to construction 
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variations; particularly bad transducers were simply not used. Since it was found that the larger transducers 
could be made with more stable characteristics, these were used exclusively on the side wall where the 
boundary layer was of sufficient thickness. 

An attempt was made to clarify a phenomenon pointed out by Willmarth^ In calibrating his trans- 
ducers, which were similar to those described here, he found that the electrical capacity between the 
ceramic and the backing could change when a pressure was applied to the sensitive surface, and that this 
could produce an undesirable component in the output of the transducer. By using a double element with both 
exposed faces grounded, he seemingly avoided this problem. Double elements could not be used here, so that 
this effect was subjected to some investigation. The tests on the transducers showed no significant capacity 
change effect, but did turn up a pyroelectric effect. A temperature variation on the surface, such as might be 
produced by a shock wave passing over the transducer, did produce an output voltage. The response time 
for this pyroelectric effect is controlled by the silver layer on the surface, and, for the elements used here, 
was about 1 millisec. The spectral region below 1 kc already had a strong contribution from the transducer 
vibrating on its mount, so that all the frequencies below 1 kc were filtered out, and this effectively elim- 
inated the temperature signal from the output. 

In order to estimate the frequency and phase response of the transducers to the imposed pressure, it 
is necessary to set up some model of the mechanical system. The pressure input can be assumed to come 
from a zero impedance source; i.e., the motion of the transducer face should not affect the pressure input. 

The experimental results related in the preceding paragraph can be interpreted as showing that the backing 
can be taken as a semi-infinite field, so that no standing waves have to be accounted for. This model 
corresponds to the classical problem of an organ pipe, forced at one end and radiating at the other. The 
model is defective inasmuch as the backing material is also forced in the neighborhood of the piezoelectric 
element by the pressure field. The analysis showed that for the material constants encountered here the 
design should give essentially uniform response, as measured by the ratio of strain to pressure for the 
frequency range in which they were used, i.e., up to 600 kc for the 0.030-in. -diameter, 0.02-in. -thick elements. 

In order to calibrate the transducers, they were located on the side wall of a shock tube and a shock 
wave was passed over them. Shock strength was determined by measuring the shock speed and temperature. 
The output signal was recorded by an oscilloscope. The response of the transducers to the shock wave was 
a sharp rise (determined mainly by the transit time of the shock over the sensitive element) and a slow 
oscillation at about 500 cps due to the mounting. A typical calibration trace, for a transducer with a 0.1-in. 
diameter and a 0.1-in. -thick PZT-5 element is shown in Fig. 5. 
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Another check on the response characteristics of the transducers was obtained by blowing a small, 
high-speed jet of air over their surfaces. Such a jet is a convenient source of an effectively continuous 
spectrum of high-frequency pressure fluctuations, and any defects in construction showed up as spikes in 
the response spectra. 

The sensitivity of the transducers ranged from 20 mv/psi to 100 mv/psi, depending on size as well 
as variations in construction. Some transducers were constructed with two or more tandem elements on the 
surface so that space -time correlations with a small spacial separation could be obtained. Insofar as could 
be determined, these elements were electrically and mechanically separate. 
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IV. ELECTRICAL EQUIPMENT 


A strain of a piezoelectric material produces a surface charge which is distributed between the 
capacity of the element and any capacity in parallel with it. The charge decays with a time constant 
determined by the total capacity and the resistance shunting the element. In order to maximize the voltage 
across the element, it is therefore essential to minimize the shunt capacity, and in order to get good low- 
frequency response, the shunt resistance should be as large as possible. 

The large resistance is obtained by using a cathode follower as the input stage of the voltage 
amplifier. The one used here was of conventional design with an input resistance of about 100 megohm. The 
cathode follower circuit included a high-pass filter with a variable low-frequency cutoff which could be set 
between 0 and 10 kc. The shunt capacity consists primarily of the lead capacity between the transducer and 
the cathode follower. It can be minimized by locating the cathode follower as close to the transducer as 
possible. An attempt to build a cathode follower that could be mounted inside the flat plate next to the 
transducer was unsuccessful because the tube microphonics in this noisy environment dominated the signal. 

It is not possible to build a solid-state circuit with sufficiently large input resistance for this purpose. The 
cathode follower was therefore located outside the tunnel, requiring about IB in. of coax cable between the 
transducer and the amplifier. 

When shaken, cables produce electrical noise, probably from the strands of the shield rubbing 
together. Commercially available low-noise Microdot cable (shield impregnated with graphite) with a capacity 
of per inch was used. The transducer capacity was normally about 15/x^f, so that a substantial 

fraction of the signal was lost in the cable. This loss can be decreased by driving the cable shield with the 
cathode follower to the voltage of the transducer; but this seemed to add more noise to the signal, and a 
passive cable was used in these tests. 

The output of the cathode follower went to a wide-band, low-noise voltage amplifier; the output of the 
amplifier went to various instruments — wave analyzers, mean-square voltmeters, etc., depending on the 
measurement to be made. A detailed description of the instrumentation is given in Ref. 7. 

The only specialized electronic instrumentation developed for these tests was a digital correlator. 

In order to measure the space — time correlation of the pressure signals on a surface, it is necessary to obtain 
the mean product of the output voltages of two transducers, where the product is formed from the output of one 
transducer at a given time and the output of the other at a fixed time interval later. To perform this operation, 
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the first technique tried was to delay one signal by a lumped circuit delay line and multiply the signals by a 
quarter-square technique. All of the correlations obtained on the flat plate were measured in this fashicn. 
This technique is limited in bandwidth since the bandwidth of a lumped circuit delay line is inversely 
proportional to the delay time increment. For large bandwidths, many small steps are required, and then the 
attenuation per step affects the output. The delay line used had 1-microsec steps and an effective bandwidth 
of about 100 kc. This bandwidth was certainly not adequate to give convincing results for the plate or for 
the sidewall boundary layer measurements. A digital correlator was developed that had a useful bandwidth 
of 200 kc and a more accurate time delay mechanism. The digital correlator was completed at about the same 
time as the sidewall mounting problem was solved so that it was used for all the sidewall measurements. 

One operation of this correlator can be described as folWs: At time zero, one signal is sampled and its 
value converted to a digit between -7 and +7 by a gate circuit. This number is stored. At time T later, the 
other signal is converted to a digit between -7 and +7. This number and the stored digit are multiplied. The 
entire operation could be repeated up to 5000 times per second, and the number of operations as well as 
the sum of the products was recorded after a convenient time interval, usually about 1 minute. The mean 
product was then computed. The delay time T could be varied continuously from 0.01 Atsec to 10 sec. 'Ihe 
correlator is described in a paper by J. Stallkamp which is to be published soon. 
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V. MEASUREMENT ERRORS 


In order to describe the statistical properties of the pressure fluctuations at the wall, two types of 
data were obtained: the temporal power spectra at a given point, and the space— time correlation between 
two points. The main sources of error in these data are the contamination of the signal with the radiated 
pressure field from the tunnel walls and the effect of the finite area of the transducer. 

Idle relative magnitude of the direct, local pressure fluctuation and the radiated pressure could be 
measured on the flat plate. When the natural transition of the boundary layer on the plate surface was down- 
stream of the transducer, a direct measure of the radiated pressure was obtained. If transition was moved 
upstream of the transducer by using the air trip but keeping the tunnel conditions fixed, the local pressure 
fluctuations were obtained. A set of these measurements is shown in Fig. 6. The radiated pressure is seen 
to be an insignificant contribution to the total signal. 

It was not possible to perform a similar experiment on the tunnel sidewall, but it is assumed that 
the same result would have been obtained. It should be noted that the two cases are not directly comparable 
since, for the flat plate, the radiation originates in the tunnel boundary layers, which are of significantly 
larger scale than the boundary layer on the plate; whereas for measurements on the sidewall, the radiating 
boundary layers and the measured boundary layer are of comparable scale. Laufer’s measurements of the 
radiation field, however, show the sound pressure level to be significantly smaller than the direct pressure 
level. 

The problem of the effect of the finite size of transducers on the measured spectra and correlation 
can be attacked experimentally by making measurements with a variety of transducers of different size and 
then extrapolating the results to zero transducer size. In order to apply this method, it is necessary to be 
able to make measurements with transducers approaching a point in size compared to the pressure scales, or 
the extrapolation will be subject to much uncertainty. An alternative approach was used here. Corcos^ has 
treated the length correction problem analytically, and has shown that with a set of assumptions that are 
consistent with the known experimental findings a usable length correction technique can be obtained. His 
results, with some modification, were used here. 

The space — time correlation function for the pressure fluctuations is defined as 
< p 2 > R ( Ax, Ay, t) = < p(x, y, t) p{x + Ax, y + Ay, t + r) > 
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where ^ ^ denotes a time average, p is the fluctuating component of the pressure at the wall, and x and y 
are Cartesian coordinates in the plane of the plate, x being aligned with the flow. It is assumed that the 
pressure field is sufficiently homogeneous in the x , y plane so that the variation of /MAx, Ay, t) with x 
and y can be neglected compared to its variation with A*, Ay and r. 

If we measure the output voltage e at a point with a transducer of sensitivity 5 volts/(sq. in.)(psi), 
the measured mean-square voltage is given by 

<e 2 > = s 2 <p 2 > //// - - it Y ~ V* 0) dxdyd^drj 

where the integration is over the surface of the transducer in both sets of variables x, y and 77. Conse- 
quently, in order to find <p 2 > from<e 2 > we need the correlation function of the pressure. The power spectrum 
of e and the space -time correlation of e are most easily accessible to measurement. Therefore, we would 
like an estimate of the integral above that is obtainable from the properties of these quantities. 

The two assumptions that give this estimate from the measured quantities are: 

1. The correlation function is a function of (x — £) and t only through the combination 

x - g - U T f where U c is a parameter, the convection speed, determined by experiment. 
This assumption permits one to obtain the x variation of the correlations (spectra) from 
the measured time spectra. When we later discuss the measured correlation functions 
of e, the justification and limitations of this assumption become apparent. 

2. The correlation function is separable in the x, y plane; i.e., R (x - y - 77, 0 ) is equal 
to Q ( x - £) G (y - 77), and, in particular, we will assume that Q (x - £) is the same 
function as G (y - 77). This assumption is made in order to get explicit results and is not 
contradicted by available experimental information. Some low-speed measurements by 
Corcos of the correlation G{y -r\) indicate that the correlation function R has essentially 
the same scale in the x and y direction. 

If we apply these assumptions to a transducer with a square sensitive area of side length 2 L, 
we obtain the result 
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<P > = 


<e 2 > 

T 2 ^ 2 



EJco) 

m 



n 2 


da) 


where co - 2tt x (frequency), E m (aj) is the measured power spectra of the voltage and is scaled such that 

- oo n 

j E (co) dco - 1, and A is the sensitive area, 4 L . Within the same framework of assumptions, the true 
0 m 

shape of the time spectrum of the pressure fluctuations at a point can be obtained by dividing the measured 
spectrurp of e at each frequency by the function 



The results for a round sensitive area are similar to these results if E (a>) falls to small values for 

m 

(ojL/U c ) < 1. With the approximation 



for {cor/U c )j rK ^ < < 1, and where /j is the Bessel function of the first kind, r is the sensitive area radius, 
and K g is the fluctuation wave number in the y direction, one obtains 
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The sensitive area of the transducer was nominally round so that the correction formula using the 
Bessel function kernel was used in reducing the data. The shape and size of the sensitive area are not 
necessarily the same as the shape and size of the ceramic element, since the mechanical working of the 
element affects the distribution of polarization. A typical example of how the correction proceeds is shown 
in Fig. 7. The measured spectrum is first divided by 


/ 


2 

1 




2 


where U is obtained from the measured correlation function. This gives a first approximation to the true 
c 

shape of the spectrum of the pressure at the measuring point. Then the level is changed by the ratio of the 
spectrum area before correction and after correction. The area under this new curve is then the correct area 
for determining <p 2 >,* <e 2> is obtained by integrating the measured spectrum <e 2 > E m (oj). 

This correction procedure gives a first approximation to the true pressure amplitude at a point but a 
somewhat rough approximation to the shape of the true time spectrum. Therefore, only those spectral shapes 
are shown in this report where the correction was not large, and for the most part they are shown uncorrected. 

As a final remark, it should be mentioned that in those few cases where sufficient data were 
available to find <p 2 > from several measurements of <e 2 > obtained with transducers with sensitive areas 
of different sizes, the correction procedure outlined above gave about the same value for <p 2 >. 
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VI. MEASUREMENTS 


A. Spectra 

The area normalized power spectra of the pressure fluctuations are shown in Fig. 8. Most of these 
data were obtained on the tunnel sidewall. The frequency coordinate used is f8/U f where 8 is the geomet- 
rical boundary layer thickness, f is the measured frequency, and U is the free-stream velocity. In this 
particular coordinate system it is apparent that^for decreasing Mach numbers, more energy becomes concentrat- 
ed at smaller nondimensional frequencies. 

The sidewall data shown are not corrected for the finite transducer size. If we assume that the data 
give the true shape of the pressure spectra, then an estimate can be obtained of the maximum transducer 
size that will yield reasonable values for <p^> without correction. The correction should be small for 
f8/U ^ < 5 on the basis of the given data. The length correction depends on the parameter 2rrfr/U c , and if 
corrections to the spectral shape should be less than 25% then 2 77 fr/ U c should be less than 1. If 

fd 277 fr r U c 1 0.6 1 

= 5 and = 1 then — = * = = 

U V 8 U 1077 1077 50 

oo c oo 

where 0.6 is a measured value of U /U . 

C oo 

Therefore, r/8 should be less than 1/50. For the sidewall boundary layer, r/8 = 1/40 at M - 4.7 and 
r/8 « 1/20 at M = 1.33. Consequently, the size effect is not particularly small for the low Mach number cases. 
Owing to the tentative nature of the length correction, however, the data are shown uncorrected. With the 
transducer sizes used, it is not possible to obtain meaningful spectral data for large frequencies, i.e., 
fS/U > 5, or get the Reynolds number effect on the spectral shapes for these frequencies. 

Some corrected flat plate measurements are included in the figure to show that the spectral shapes 
are roughly the same as for the wall measurements. 

If we assume that the pressure field is convected past the transducer as a frozen pattern with a 
velocity U (the evidence for this assumption is covered in the following section) then it is possible to 
relate the measured frequency to the spacial scales of the fluctuation field. In particular, the integral scale 
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of the pressure correlation R(Ax, 0, 0), defined as 


L = f R(kx, 0, 0) d(Ax) 

P 


is simply related to the intercept of the pressure spectrum on the /= 0 axis. For the coordinate used, i.e. 
fS/V oo , the result is 


4 1 / 


Up o 


E_ (0) 


The values for this ratio range from 0.16 at M = 1.33 to 0.057 at M — 4.54. Since this variation is approximately 
the same as the variation in the ratio of the boundary layer momentum thickness 6 to the geometrical thickness, 
the spectra are plotted in Fig. 9 using the frequency coordinate fO/U c - This plotting method brings all the 
data except those at M = 1.33 into approximate coincidence. Since the effect of the finite transducer size is 
to increase the apparent scale, the M = 1.33 data, where this error is the largest, might be showing this 
effect. However, the M = 2.00 data would have about the same length correction and the scale is significantly 
smaller than at M =1.33. The conclusion is, therefore, that the relative scale does increase for Mach 
numbers less than 2. 

A measurement by Morkovin (9 * of the spectrum of the velocity fluctuations in a Mach number 1.76 
boundary layer shows that the length scale of the velocity fluctuations here is about the same as for an 
incompressible boundary layer, i.e., L /S = 0.4. Therefore, even at M = 1.33 the pressure field has a 
significantly smaller length scale than the velocity field. 

Figure 10 shows measurements of hot wire spectra obtained in the sidewall boundary layers. The 
wire operating conditions were such that the primary sensitivity was to temperature fluctuations. In agreement 
with Morkovin’s results.the length scales of the temperature fluctuations are 0.1 of the boundary layer thick- 
ness. The average pressure spectrum for M > 2 is sketched in the figure to show how the pressure spectra 
compare with the temperature spectra. Also shown is the velocity spectrum obtained by Klebanoff* 1 in an 
incompressible boundary layer. The results of Morkovin indicate that the compressible velocity spectra are 
about the same as the incompressible spectrum measured by Klebanoff. 
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B. Correlations 


The space — time correlations considered here are 

<p > R(Ax, r) = <p(x, y, t) p(x + Ax, y, t + r)> 


where x is the stream-line direction. As is known from subsonic results (Willmarth), the shape of the 
correlation curve obtained for a fixed longitudinal spacing of the transducers and a variable time delay 
between the two signals displays a maximum for a particular value of the time delay, T max • A velocity, the 
so-called convection velocity, can be constructed by dividing the spacing Ax by the time T and this 
velocity can be interpreted as the speed at which the pressure field is moving along the wall. The correlation 
function can be characterized by this speed and by the rate at which the value of R{£, t) decreases for 
increased spacing of the transducers. Only the first variation is included in the length correction calculations. 
Typical measured correlation functions are shown in Fig. 11 and 12. The results in Fig. 11 were obtained 
on the sidewall with the electronic correlator; the results in Fig. 12 were obtained on the flat plate using 
the lumped circuit delay line. 

The convection speed Ax/t obtained from these data is shown in Fig. 13 as a function of Mach 
number. The speeds shown were computed only from sidewall data, obtained with the electronic correlator. 
Willmarth’s subsonic data are also shown. At M = 2.00, the tunnel Reynolds number was changed by a factor 
of 10, and it was found that there was no measurable change in the correlation shape or in the convection 
speed over this range of Reynolds number. 

The rate of decay of the peak values of the correlation coefficient for various Mach numbers is 
illustrated in Fig. 14a, where the peak value is plotted vs the nondimensional distance Ax/3. The absolute 
value of the correlation was measured by extending the time delay in the measurement until the correlation 
had dropped to a constant value, which in the absence of extraneous signals should have been zero. The 
peak value was assumed to be the difference between this level value and the measured peak. 

The data shown for Ax/5 less than 0.6 were obtained on the sidewall and for Ax/5 > 0.6 on the 
flat plate. Assuming that these data are comparable, curves of the maximum value as a function of distance 
are sketched in. The general shape of the curve is a sharp drop from the value 1 at Ax = 0 to about the 
value of 1/2 at Ax/5 ~ 0.2, and then a more gentle variation for larger Ax. There is no outstanding effect of 
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Mach number on the decay of the peak value, but careful examination of the data seems to indicate a trend for 

jJe 

the peak to fall more rapidly for smaller Jf M . The peak values are plotted against Ax/S in Fig. 14b, which 
gives a little better collapse of the data for large M x and Ax. 

C. Fluctuation Levels 

The determination of the fluctuation levels from the calibrated transducers requires the integration of 
the measured power spectra. This method is more desirable than the direct measurement of the mean-square 
voltage because energy associated with the transducer can be removed from the spectra, and a length 
correction can be applied to the spectra before integration to obtain a more accurate result. 

In the equilibrium flat-plate boundary layers considered here, it is to be expected that the rms 
pressure fluctuation level divided by some characteristic pressure, such as the local free-stream dynamic 
pressure, should be a function only of the Mach number and a Reynolds number based on the local boundary 
layer thickness. The Reynolds number variation of <p 2 > was investigated primarily on the flat plate, since 
the Reynolds number could be varied over a wide range there by varying not only the tunnel conditions but 
also the location of the transducer on the plate. 

When the mean pressure level is varied and the spectra are measured at one point, the shape and 
magnitude of the measured spectra change for two reasons. As the mean pressure is increased, the boundary 
layer becomes thinner, and consequently the energy associated with the large eddies spreads to higher 
frequencies. At the same time, the spectra are affected by the change in Reynolds number. The spread of the 
energy to higher frequencies changes the length correction for different Reynolds numbers, and in order to 
obtain the pure Reynolds number effect, the length correction must be introduced. 

The data obtained from the flat plate were corrected, as explained earlier. The size of the correction 
on the measured <e 2 > was as large as a factor of 2 in some cases. The corrected data are shown in Fig. 15. 
Owing to the size of the length correction, the relative values of <p 2 > obtained from different spectra at the 
same Mach number but at different Reynolds numbers are more accurate than the absolute values. Therefore, 
only the values of p'/p' (Re ^ = 10*) are shown in Fig. 15. The absolute level was determined from the side- 
wall measurements, where the corrections were much smaller. It is seen that for all Mach numbers the 
corrected levels of 'p/q have a variation with Reynolds number close to p/q ~ Re^’^ 9 which suggests that 
if the skin friction were used as the characteristic pressure, the effect of Reynolds number on p would be 
accounted for. 
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The pressure fluctuation levels obtained on the sidewall are shown in Fig. 16, where r , the mean 
skin friction at the wall, is used to remove the Reynolds number influence. The magnitude of p '/t changes 
only little with Mach number above > 2; was computed using the compilation of Van Driest^ 

The subsonic results of Willmarth and Harrison are plotted on the same scale and are seen to be 
significantly lower than the supersonic results. In order to see if our technique gave the same results at low 
speeds, one measurement was made on the flat plate with the tunnel choked near the diffuser so that there 
was subsonic flow over the flat plate. The sharp leading edge had to be modified with a wood filler to 
prevent the flow from separating. One pressure fluctuation measurement, also plotted on the curve, was 
obtained at M = 0.6. This point is higher than Willmarth’s data, but does indicate a trend of smaller pY T 
at lower speeds when compared to the supersonic data. 
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VII. DISCUSSION 


The data that exist concerning the structure of the turbulence in supersonic boundary layers show 
that the velocity field is not strongly influenced by Mach number effects, at least in the outer regions of the 
boundary layer^ The measurements reported here of the wall pressure field of a supersonic boundary layer 
show that the root-mean-square level of the pressure divided by the local shear is higher than the subsonic 
value but remains reasonably constant above Mach number 2. The scale of the pressure fluctuations, however, 
changes in relation to the scale of the velocity field. The length scale for the pressure field above Mach 
number 2 is proportional to the momentum thickness of the boundary layer, but it is not clear what physical 
significance can be ascribed to this fact. The momentum thickness is a defined quantity whose gradient is 
related to the local skin friction but whose numerical value is not meaningful with respect to the dynamics 
of the turbulence. 

The most probable explanation for the change in the scale of the pressure field in contrast to the 
velocity field is that the existence of regions of relatively supersonic mean flow affects the structure of the 
pressure field within the boundary layer. A turbulent boundary layer consists of a relatively disordered 
collection of vortices embedded in a region of mean vorticity. A turbulent eddy consists of a rotational part 
and an associated quasi-potential velocity field. The rotational core is swept along at the local velocity and 
the potential field follows the core. This potential part extends in all directions, and, in particular, can extend 
into regions of the flow where the mean velocity is supersonic with respect to the core velocity. The presence 
of these relatively supersonic regions could attenuate the effect on the wall of a large eddy in the outer 
region of the boundary layer and could amplify the effect on the wall of an eddy near the wall. Such a physical 
process would decrease the effective length scale of the pressure fluctuations on the wall as well as decrease 
the convection velocity by limiting the range of mean velocities of the eddies that produce the major effect at 
the wall. The fluctuating velocity field would not be strongly influenced by such a mechanism since a large 
part of the measured velocity fluctuations arise from the displacement of the mean velocity gradient and this 
part of the velocity fluctuation is not coupled to the pressure field. 

The maximum mean velocity of an eddy whose influence is subsonic at the wall is computed in an 
appendix. An adiabatic boundary layer is assumed. If this velocity is converted to distance from the wall 
assuming a 1/7 power velocity profile, z/h for such an eddy varies from 0.39 at M = 1.33 to 0.06 at M = 2.00. 
Therefore, the range of mean velocities where an eddy has a subsonic influence on the wall decreases rapidly 
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as the Mach number increases, and already at M ~ 2.00 the corresponding region of flow is only a small 
fraction of the boundary layer thickness. 

The space — time correlation along the stream direction changes rapidly in magnitude; in fact, it has 
fallen to half its initial peak value within 1/5 of a boundary layer thickness. This rapid falloff indicates that 
the eddies contributing to the pressure at the wall have widely differing velocities or that the eddies them- 
selves vary rapidly in a coordinate system fixed to their centers. Since measurements of Favre^ 2 ^ of the 
space— time correlation of the velocity fluctuations in a subsonic boundary layer have shown that the tur- 
bulence pattern — at least away from the wall — moves with the local velocity and does not vary rapidly in 
distances comparable to the boundary layer thickness, it is probable that the pressure correlations on the 
wall — and perhaps also the velocity correlations near the wall — decrease to small values within a short 
distance because of the wide range of mean velocities of the large eddies that add up to give the wall effects. 

A useful theory for the pressure fluctuations would have to predict not only the fluctuation levels 
but also the properties of the correlation function. If the ideas expressed above are correct, this would imply 
that some knowledge of the relationship of the large eddy motion in the turbulence to the pressure field at 
the wall would be required, as well as a knowledge of the essential properties of the turbulent eddies. For 
example, what is the distribution of the location of the centers of the eddies, i.e., that part of the organized 
eddy motion that moves with the local mean speed? The currently available theories bypass such problems 
with ad hoc assumptions, and important information is lost. 

Most of the current information on the properties of turbulent boundary layers has been obtained with 
a hot-wire anemometer. Therefore, it is of interest to examine the way in which the results of this work can 
affect the interpretation of a hot-wire signal obtained from a supersonic boundary layer. It has been shown^ 1 ^ 
that the output of a hot wire in a supersonic boundary layer can be decomposed into a contribution by a 
velocity and by a temperature or density fluctuation if one can assume that the pressure fluctuation level is 
much smaller than the levels of the other fluctuation; i.e., 

P P T U 

P P T U 

Since the experiments reported here show for M > 2, p/r^ = 5 at the wall, we can write 
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Therefore, 

p/p « 2.5 y A/ 2 Cy 


A typical value for C y is 10 ^ and y = 1.4, so that 

~ ~ 3.5 x 10" 3 A/ 2 

P 

For Af^ = 5, p /p = 8.75 x 10~ 2 = 8.75%. This value is comparable to the velocity fluctuation level in a 
boundary layer, so it appears that for M in the neighborhood of 5 or larger the simple constant pressure 
interpretation of the data is no longer possible. Any results obtained for these large Mach numbers will have 
to incorporate some knowledge of the pressure field and its relation to the turbulent velocity field. 
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VIII. CONCLUSION 


Measurements of the wall pressure field associated with a supersonic turbulent boundary layer have 
been obtained for flow Mach numbers up to 5.0. The major effect of increasing the Mach number is to decrease 
the length scale of the pressure field. The integral scale of the wall pressure fluctuations changes from 16% 
of the boundary layer thickness at M = 1.33 to 0.06% at M = 4.54. The rms level of the pressures is propor- 
tional to the local mean shear for all Mach numbers,and the proportionality constant changes from about 3 for 
subsonic boundary layers to about 5 for Mach numbers greater than 2. 

The convection speed characterizing the pressure correlation on the wall decreases with increasing 
Mach number, and within measurement accuracy is independent of Reynolds number. The maximum value of 
the pressure correlation falls off rapidly with increased spacial separation of the measuring points, and the 
peak is down to a value of 1/2 at about A* = 0.2 8. The rapid decrease of the correlation implies that the 
eddies contributing to the wall pressure move with a large range in speeds, so that the convection speed 
represents some weighted average of these speeds. 

The main limitation to the data reported here lies in the uncertain nature of the length correction, 
but this limit can only be removed if a better grasp of the structure of the pressure field throughout the 
boundary layer is obtained. 
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NOMENCLATURE 

a sound velocity 

A area of sensitive region of pressure pickup 
C y local skin friction coefficient = T^/ q 

d diameter of round transducer 

e fluctuating component of transducer output voltage 
E m measured power spectra of e, normalized so that J E m (co) d to = 1 
E normalized power spectra of hot-wire output 

£’ 1 E divided by J %(cor/U )/(cor/U 

i m j i c c 

f frequency, cps 

&2 wave number in y direction (1/in.) 

L one-half the length of the side of a square transducer 
Lp integral scale of the wall pressure field 

M Mach number of free stream 

oa 

p fluctuating component of pressure on surface, p mean pressure 
<p 2 > mean square value of p 

'p root-mean-square value of p = V<p 2 > 

2 

q dynamic pressure in free stream = 1/2 U ^ 

r radius of a round transducer 

Re Reynolds number, Re ^ = U ^S/v^, Re s * = 8 /v^, Re & ^ oc p! v 00 

R peak value of R{r) for a particular Ax 

max r 

/?(Ax, Ay, t) space — time correlation coefficient of pressure fluctuations on surface 
<p 2 > R{ Ax, Ay, r) = <p(x, y, t) p(x + Ax, y + Ay, t + r)> 

R(r) space — time correlation coefficient for output voltage of two transducers separated by a 
distance Ax-<e 2 > R(r) = <e(x, t) &{x + Ax, t + t)> 

s sensitivity of pressure transducer per unit area (volts/psi(inch) 2 ) 

_ 

T, T mean static temperature and rms value of fluctuations from mean 
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NOMENCLATURE (Cont’d) 


v . 

'V. 

u, u 


u 


U 0> U l 


Ax 

y 

Ay 

Z 

y 

8 


v 

e 

€ 

p 

p 

T 


r 

max 


r 

w 


co 


< > 


convection speed = Ax/i; where T is at location of R 

1 max 

mean velocity and rms value of fluctuations from mean 
free-stream velocity 

mean velocities in boundary layer such that — Uq = a^ 9 where is sound speed at 
point where velocity is U ^ 

Cartesian coordinate in plane of solid surface aligned with mean stream direction 
separation of two transducers in x direction 

Cartesian coordinate in plane of solid surface perpendicular to x 
separation of two transducers in y direction 
Cartesian coordinate perpendicular to solid surface 
ratio of specific heats of a gas, here taken at 1.4 for air 

thickness of boundary layer, i.e., distance from solid surface where velocity attains 
free-stream value 

displacement thickness of boundary layer 

alternative notation for y, i.e., another coordinate in y direction 

momentum thickness of boundary layer 

a coordinate in x direction 

mean density of air 

rms value of fluctuations from p 

time increment 

value of T for which /?(t) has maximum value 
mean skin friction at wall, psi 
angular frequency, co - 2rr( 
time average 
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APPENDIX 

Maximum Speed of an Eddy with Subsonic Influence at the Wall 

The data presented in this report make it clear that the convection speed represents not some unique 
velocity for an eddy pattern, but instead is the weighted average over many eddy velocities. Still, it is of 
some interest to examine the relationship of this speed to the various sonic velocities of importance in a 
boundary layer. Here we refer to the sonic lines relative to an eddy moving with some velocity i f 0 with 
respect to the wall, i.e., the boundary between the region where the eddy has a subsonic influence and the 
region where the eddy has a supersonic influence. Such a sonic line satisfies - Uq = flp where U j is the 
local mean velocity on the sonic line. If we assume that we are dealing with a perfect gas flow with constant 
specific heats, then for an adiabatic boundary layer the energy equation is C p T + (U 2 / 2) = C p T Q (C p is the 
specific heat at constant pressure). It follows that the condition relating U l and U Q is 



V 

oo 



4 U 0 U 1 2 / ^0 \ 2 7-1 2 1 

7+1 f /2 7 + 1 \ U" / 7+1 7 + 1^/2 

OO OO 


where U and M are the free-stream velocity and Mach number, respectively. 

OO oo 

Of interest is the eddy velocity Uq where the wall is a sonic line; i.e., U ^/U ^ = 0, or the maximum 
speed at which the eddy has a subsonic influence at the wall. We obtain 



2 


0.2 


M* 


if y = 1.4. For an eddy moving with this velocity, the region of subsonic influence extends away from the 
wall to a region of velocity 


u i 


u 


oo 


2.78 
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These two curves are plotted in Fig. A-l, along with the measured values of the convection speed. Of 
interest is the fact that even at M = <*>, under the condition postulated, an eddy with a subsonic influence at 
the wall has a subsonic influence in the boundary layer out to the point where U^/U^ = 0.75, a finite fraction 
of the boundary layer thickness. This serves to point up the fact that some care has to be taken in making 
assumptions with respect to the limit of M -► oo. 

Phillips^^ assumed that the width of the subsonic region, in an analogous problem, shrinks to 
zero essentially as 1/M as M -»■ oo. This implies a similarity of the Mach number field as M -» oo, which is not 
true if Tq = constant is a constraint. 

The experimental values for the convection speed are approximately equal to the eddy velocity that 
has a sonic influence on the wall up to M = 3. Above M = 3, the convection speed levels off at U (; /U oo = 0.6, 
and the eddies primarily responsible for the wall pressure move supersonically with respect to the wall. 
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PRESSURE TRANSDUCERS 


Fig. 1. Transducer installation on flat plate 
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Fig. 2. Transducer installation on tunnel side wall 
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Fig. 4. Photograph of transducers installed 
in mounting disk 


Fig. 3. Drawing of pressure transducer 


Fig. 5. Shock tube response of transducer 
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Fig. 10. Spectra of a hot-wire signal near the center of the boundary layer 
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Fig. 12. Measured space— time correlations 
on the flat plate 



Fig. 13. The convection speed ratio 
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PEAK VALUE OF CORRELATION FOR VARIOUS Ax 
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Fig. 14. Peak values of the correlation for 
various transducer spacings 
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